Notes on anelastic effects and thermal noise in 
suspensions of test masses in interferometric 
gravitational- wave detectors 

Constantin Brif 
LIGO Project, California Institute of Technology, Pasadena, CA 91125 



Contents 



1 Introduction 



2 Fluctuation-dissipation theoremi 



2.1 Example: Damped harmonic oscillator 



3 Anelasticity of solids 



3.1 The complex Young modulus and the loss function 



3.2 Simple models 

|3^.1 Perfect elastic solid 



J.2.2 Maxwell sohd 



j.2.3 Voigt-Kelvin solid 



|3.2.4 Standard anelastic solid 



3.3 Boltzmann's superposition principle 



3.3.1 Example: Standard anelastic solid 



4 Calculation of the thermal noise spectrum for a pendulum suspension 



2 

3 

3 

3 
4 
4 
4 
5 
5 
6 



8 



4.1 The direct approach 



4.2 The normal- mode decomposition 



9 

4.3 Modes of a pendulum suspension 10 

4.3.1 The pendulum mode 10 

12 



4.3.2 The violin modes 



5 Experiments on anelasticity effects for pendulum suspension^ 



13 

5.1 Basic types of experiments 13 

5^2 Internal losses in wire materials! 14 



6 Conclusions 



Appendix: Correlation function and spectral density 



Bibliography 



15 
15 
16 



1 Introduction 



The thermal noise is expected to be one of the main hmiting factors on the sensitivity of interfero- 
metric gravitational- wave detectors like LIGO and VIRGO. Thermal fluctuations of internal modes 
of the interferometer's test masses and of suspension modes will dominate the noise spectrum at the 
important frequency range between 50 and 200 Hz (seismic noise and photon shot noise dominate for 
lower and higher frequencies, respectively). It is important to note that off-resonant thermal noise 
level in high-quality systems is so low that it is unobservable in table-top experiments. Therefore, 
predictions of the thermal-noise spectrum in LIGO are based on a combination of theoretical models 
(with the fluctuation-dissipation theorem of statistical mechanics serving as a basis) and experimental 
measurements of quality factors of systems and materials involved. It is assumed that losses in the test 
masses and suspensions will occur mainly due to internal friction in their materials, which is related 
to anelasticity effects in solids. 

These informal notes comprise some basic results on the theory of anelasticity and thermal noise 
in pendulum suspensions. This collection is by no means complete and focus on aspects which are 
of interest for the author. The original results can be found in a number of books, research papers, 
and theses. Some of these sources are listed in a short bibliography at the end of the present text; a 
list of research papers (since 1990) devoted to various aspects of the thermal noise in interferometric 
gravitational- wave detectors was prepared by the author and is available as a separate document. 



2 Fluctuation-dissipation theorem 

Consider a linear one-dimensional mechanical system with coordinate x{t). If a force F{t) acts on the 
system, than in the frequency domain the force and the coordinate are related by 

x{io) = H{uj)F{uj), (2.1) 

where H{uj) is the system response function. Then the spectral densities (see Appendix) are related 

by 

S^ioj) = \H{oj)\^Sf{oj). (2.2) 

The impedance of the system is defined as Z{uj) = F{u>)/v{u)) = F{Lj)/[iLux{u))]. Therefore, the 
impedance and the response function are related by Z{uj) = l/[iujH{uj)]. 

If the system is in equilibrium with the thermal bath of temperature T, then the fluctuation- 
dissipation theorem (FDT) says that the spectral density of the thermal force is 

Sf{Lo)=4kBTReZ{uj), (2.3) 

where ks is the Boltzmann constant. The form ( |2.3D of the FDT is valid in the classical regime, when 
the thermal energy ksT is much larger than the energy quantum hui. Using the FDT, one readily 
obtains the thermal noise spectrum 

S'J\co) = l^Rey(a;), (2.4) 

where Y{lo) = 1/Z{uj) is the admittance and IieY{uj) = KeZ{uj)/\Z{uj)\'^ is the conductance. The 
FDT is the basis for calculations of the thermal noise spectrum in interferometric gravitational-wave 
detectors. 
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2.1 Example: Damped hcirmonic oscillator 

Consider a damping harmonic oscillator of mass m, spring constant k, and damping constant 7. The 
equation of motion is 

mx + + kx = F(t). (2-5) 
In the frequency domain this can be written as 

{-mw^ + i'jLO + k)x{u) = F{lo). (2.6) 

The impedance of this system is Z{(jj) = ^ + i{mui — k/u)). Then the FDT gives the spectrum densities: 

Sf{u) = 4kBT^, (a;) = -. /^fl^ , 2 - (2-7) 

3 Anelasticity of solids 

The FDT means that if a system has no dissipation channel, thermal fluctuations will be zero. For an 
ideal elastic spring without friction, ReZ(a;) = 0, and there are no fluctuations: -S'*'^(a;) = 0. Deviations 
of solids from the ideal elastic behavior (anelasticity) will result in internal friction (dissipation) and 
related thermal noise. For gravitational- wave detectors like LIGO, the test masses will be highly 
isolated, so the internal friction in materials of which the masses and their suspensions are made is 
believed to be the main source of dissipation and thermal noise. 

3.1 The complex Young modulus and the loss function 

Deformations of solids are usually described in terms of stress a and strain e (equivalents of mechanical 
restoring spring force Fg and displacement x, respectively). Perfect elastic solids satisfy Hooke's law 

a{t) = Ee{t), (3.1) 

where E is the (constant) Young modulus (an equivalent of the spring constant k). Anelasticity can be 
described by introducing the complex Young modulus (or the complex spring constant in a mechanical 
model). This is done in the frequency domain: 

E(.) = ^, *(u.) = (3.2) 

If an external force F{t) acts on a point mass m attached to such an anelastic spring, then the equation 
of motion in the frequency domain is 

[-ma;2 + k{ijj)]x{ijo) = F{uj). (3.3) 

The impedance of this system is 

and ReZ{u) = (l/a;)Im A;(a;). Now, the FDT theorem gives the thermal noise spectrum: 

^'^^'^'^ = kRLO (l-ma;7fci)2 + ^2- (3-5) 
3 



Here, kR{uj) = Rek{Lo), and 

^(-) = (3-6) 

is the so-called loss function. Note that cf) = tan (5, where S is the angle by which strain lags behind 
stress. The loss function is a measure of the energy dissipation in the system. The rate at which 
energy is dissipated is FgX. Then the energy dissipated per cycle by an anelastic spring is 



AS = {2'k/uj)F,x. (3.7) 
Taking Fg = Fq cos tot and x = xq cos(c<;t — 6), one finds 

AS = TrxoFosm6. (3.8) 
If 5 is small than the total energy of spring vibration is S = ^xoFq. Then for 6 <^1, one obtains 

For small (p (which is usually the case for the internal friction in materials used in detectors like 
LIGO), it is customary to neglect the frequency dependence of k^. Then one can write k{u) = 
k[l + i4>{uj)], where k = itiloq is a constant (and luq is the resonant frequency). Though this is a good 
approximation for many practical reasons, in general kji must be frequency-dependent because real 
and imaginary parts of k{ui) are related via the Kramers-Kronig relations. 

3.2 Simple models 

Here we consider some simple models of anelasticity in solids. Neither of them gives a full description 
of the behavior of a real material, but nevertheless they are useful from the didactic point of view. 

3.2.1 Perfect elastic solid 

The mechanical model of perfect clastic solid is a lossless spring. In this case a = Ee, so cj) = 0. There 
is no dissipation and no thermal noise. 

3.2.2 Maxwell solid 

The mechanical model of Maxwell solid is a lossless spring in series with a dashpot. The dashpot 
provides a source of viscous friction with a = rje. Then for Maxwell solid stress and strain are related 
by equation 

e = E-^a + rj-^a. (3.10) 

This equation shows that for a constant strain, stress decays exponentially. On the other hand, for a 
constant stress, strain increases linearly, which is a very wrong description for crystalline solids. Going 
to the frequency domain, one obtains 

A _ ^(^) _ ^'^^^ _ ^^^^-^ + i^rjE^ 
^'^^ " e(a;) ~ E + iurj ~ E^ + cjW ^ ^ ^ 

and (piio) = E/{r]Lo). 
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3.2.3 Voigt-Kelvin solid 

The mechanical model corresponding to Voigt-Kelvin anelastic solid consists of a lossless spring and a 
dashpot in parallel, which corresponds to a damped harmonic oscillator. The relation between stress 
and strain reads 

r]e + Ee = a. (3.12) 

For a constant stress (Tq, strain changes exponentially with the decay time r]/E from its initial value eg 
to the equilibrium (Hooke) value ctq/E. For a constant strain, stress is also constant, like in Hooke's 
law. This is a good description for materials like cork, but it is not suitable for metals. In the frequency 
domain, one has 

E{uj) =E + irjLO, (l){oj) = (j]/E)uj. (3.13) 



Substituting this (j) into Eq. (^), we find 

(mw^ — Ey + rj'^uj'^ 

This is the same as Eq. (p.TD for a damped harmonic oscillator with E ^ k and rj ^ ^. 



3.2.4 Standard anelastic solid 

Though the model of standard anelastic solid (SAS) does not gives a complete account of properties 
of real metals, it describes quite well basic mechanisms responsible for anelastic effects. In fact, if a 
dissipation mechanism has characteristic relaxation times for strain upon a constant stress and for 
stress upon a constant strain, then the SAS model gives an adequate description. The corresponding 
mechanical model consists of a spring in parallel with a Maxwell unit (which is a spring in series with a 
dashpot). Let Ei and E2 be the Young moduli of the separate spring and of the spring in the Maxwell 
unit, respectively, and rj be the dashpot constant, as usual. Then stress and strain are related by the 
following equation: 

^a + a ='^^e + {Ei + E2)e. (3.15) 

For a constant strain eo, stress decays exponentially from its initial value ctq to the equilibrium (Hooke) 
value EieQ-. 

a{t) = Eieo + {ao - Eieo)e-'/^% T, = r^/E2. (3.16) 

Analogously, for a constant stress ctq; strain decays exponentially from its initial value eo to the 
equilibrium (Hooke) value (Jq/Ei: 



•^0 , f, ^o\„-t/r<^ _ E1 + E2 



Then Eq. ( 3.15| ) can be rewritten in the following form 

a + T,& = ER{e + T^e), (3.18) 

where Efi = Ei is called the relaxed Young modulus. Transforming to the frequency domain, one 
obtains 

(1 + iiOT^)a{uj) = Eji{l + iLOT^)e{uj). (3.19) 
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Then the complex Young modulus is given by 



Ei.) = En\±^ = ^^ (l+^^Yj + M^^-^^) . (3.20) 

1 + ZWTg 1 + U^'t^ 



It is easy to see that 




E{u)^{ '^'^l (3.21) 



= = Ar^, (3.22) 



where = Ei + E2 is called the unrelaxed Young modulus. The loss function has the form 

1 + LO^TaTe 1 + W^f 2 

where 

T = VTaTe = Te\—, A= ^.^^ = -. (3.23) 

One sees that (p oa lo for tof <^ 1 and (p oc for u;f 3> 1. The loss function has its maximum 
^ = A/2ata; = f~^. This is called the Debye peak. This behavior is characteristic for processes 
with exponential relaxation of stress and strain, f is the characteristic relaxation time and A is the 
relaxation strength. 

Thermoelastic damping mechanism Zener pointed out that the SAS model with 



is suitable for describing processes in which the relaxation of stress and strain is related to a diffusion 
process. One example of such a process is the so-called thermoelastic damping. Consider a specimen 
which is a subject to a deformation in such a way that one part of it expands and the other contracts 
(e.g., a wire of a pendulum which bends near the top while the pendulum swings). The temperature 
increases in the contracted part and decreases in the expanded part. The resulting thermal diffusion 
leads to the dissipation of energy. This anelastic effect can be described by the SAS model with the 
thermal relaxation strength and relaxation time given by 

EuTa^ _ £ 

where T is the temperature, a is the linear thermal expansion coefHcient, Cy is the specific heat per 
unit volume, d is the characteristic distance heat must flow, and D is the thermal diffusion coefficient, 
D = g/Cy, where g is the thermal conductivity. For a cylindrical wire of diameter d, the frequency of 
the Debye peak is 

/=^-2.6|. (3.20) 
3.3 Boltzmann's superposition principle 

While the SAS has certain general features in common with actual solids, it does not reproduce 
precisely the behavior of any real metal. Simple models considered above can be generalized by a 
theory which only assumes that the relation between stress and strain is linear. This assumption 
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was expressed by Boltzmann in the form of a superposition principle: If the deformation xi{t) was 
produced by the force Fi{t) and the deformation X2{t) was produced by the force F2{t), then the force 
Fi(t) + F2{t) will produce the deformation xi(t) + X2{t). On the other hand, the deformation can be 
regarded as the independent variable. In this case the superposition principle states: If the force Fi (t) 
is required to produce the deformation xi{t) and the force F2{t) is required to produce the deformation 
X2{t), then the force Fi{t) + F2{t) will be required to produce the deformation xi{t) + X2{t). 

Let us introduce the quantity A(t) which is called the creep function and is the deformation resulting 
from the sudden application at t = of a constant force of magnitude unity. During an infinitesimal 
interval from t to t + dt, the applied force F{t) can be approximated by a constant force of magnitude 
Fdt. Then the superposition principle gives 

x{t) = f \{t - t')F{t')dt'. (3.27) 



Conversely, we may regard the deformation as a specified function of time. Let us define the quantity 
K{t) which is called the stress function and is the force which must be applied in order to produce 
the step-function deformation x{t) = @{t) (here &{t) is 1 for t > and for t < 0). Then the linear 
relationship is 

F{t) = [ K{t- t')±{t')dt'. (3.28) 



The relation between the creep function and the strain function is rather complicated; in general they 
satisfy the following inequality 

\{t)K{t) < I. (3.29) 

For constant K{t) = k we recover Hooke's law F{t) = kx{t) and then A(t) = k^^. Integrating by parts 
in Eq. (3.28), we obtain another expression of the superposition principle, 



F{t) = f f{t - t')x{t')dt', (3.30) 

where 

f{t) = k{0) 5 {t) + k{t). (3.31) 

The relationship between the force and the deformation becomes very simple in the frequency 
domain. Toward this end we introduce the functions 

fp{t) = f{t)e{t), Kp(t) = K(t)e(t), Xp{t) = x{t)@{t), (3.32) 

which are zero for t < 0. Using these functions, one can expand the upper integration limit in 
Eqs. ( pTl) , ( 13^2^ ), and (|33o| ) to oo. Then we just can use the fact that a convolution in the time 



domain is a product in the frequency domain. This gives 

F{uj) = iujKp{uj)x{ijj) = fp{uj)x{uj). (3.33) 

Thus the Fourier transform of the stress function is simply related to the complex spring constant of 
Eq. (P): 

k{u) = fp{ijj) = iujKp{ijj). (3.34) 
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3.3.1 Example: Standard anelastic solid 

For the SAS the stress function is given by 

K{t) = Er + {Eu - ER)e-'/^^ . (3.35) 

It is straightforward to see that this function leads to the first-order differential equation of the form 
( |3.18D . Then we find the function f{t), 

fit) = Eu5{t) - En^lR^-t/r. ^ (3.36) 

and the complex string constant, 



oo 



k{u^)=l f{t)e-'^'dt = Eu-^^^—^. (3.37) 



This can be rewritten in the form 



k{.) = En'-^ = E^ i^+-'rr^)+Mr.-r,) ^ ^^^^^^ 
1 + luiTf 1 + w^r/ 



which coincides with Eq. ( ^.20 ). 



4 Calculation of the thermal noise spectrum for a pendulum sus- 
pension 

For a point mass m attached to an anelastic spring with the complex spring constant k{uj), we found 
a simple result 

k{ijj) — vnuP' 



Z(a;) 



iijj 



which can be used in the FDT to derive the thermal noise spectrum S'*'^(u;) as given by Eq. (3.5). 
However, the question is how to find the thermal noise spectrum for more complicated systems, e.g., 
for pendulum suspensions of test masses in interferometric gravitational- wave detectors like LIGO. 

In the literature we can find two different approaches: the "direct" application of the FDT to 
the whole system and the method which is based on decomposing a complicated system into a set of 
normal modes. Below, we describe briefly both of these approaches. 

4.1 The direct approach 

In brief, this method can be described as follows. First, one should write equations of motion for the 
whole system and find the impedance ZiuS). Then the FDT provides the thermal noise spectrum: 

5th(^) = l^Re[l/Z(..)]. (4.1) 

The impedance Z{ui) contains the information about resonances of the system. The dissipation enters 
by taking the Young moduli of the materials to be complex: 

E{uj) = [1 + Re E{uj), (4.2) 

or, for simplicity, i?(a;) = Eq\V + ii?i>(a;)], where £"0 is a constant. The loss function ^{10) is obtained 
from experiments on the anelasticity of materials used in the system (e.g., on the suspension wires). 
Of course, the resulting noise spectrum 5*^(u;) depends very much on what form of 4>{oS) is used. 
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4.2 The normal-mode decomposition 

The normal-mode decomposition is a more traditional approach. Consider, for example, an one- 
dimensional system of linear mass density p{z) and total length L, which is described in terms of the 
normal modes ipn{z)- These modes satisfy the orthonormality relation, 

L 

p{z)lpm{z)lp„{z)dz = dmn, (4.3) 



and an arbitrary displacement x{z,t) can be decomposed as 

x{z,t) = Y,Mz)Qn{t). (4.4) 

n 

Here, qn{t) are the mode coordinates which satisfy 

qn + ujlqn = Fn{t), (4.5) 
where w„ are the resonance frequencies of the modes, and 

Fn{t) = ["^ f{z,t)Mz)dz (4.6) 
Jo 

is the generalized force produced by the force density f{z,t) applied to the system. 

This decomposition effectively replaces the complicated system by a collection of oscillators, and 
each of them satisfies 

[-u;' + u;l{uj)]qn{u;) = Fn{uj). (4.7) 
The dissipation is included by taking 

ivliu;) = ojI[1 + iMiv)], (4.8) 
where (j)n{^) are the loss functions. Then we can write 

Qniu^) = . , (4.9) 

Let us assume that the force is applied at the end of the system z = L, such that f{z,t) = 
F{t)5{z — L). Then the generalized forces are Fn{t) = F(t)'ipn{L), and we can substitute Eq. ( [4.9| ) into 
the Fourier transform of Eq. ( [4.4| ) to obtain 

x{L,u)=x{u;)=T 5 ^$^^^^-^-—F{lo). (4.10) 

This gives the admittance of the system in the form 

Then the FDT can be used to obtain the spectral density of thermal fluctuations at z = L: 
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This can be written as a sum 



over the contributions 



5fH = ^5^na;) (4.13) 



> u {l-mnuVkn? + <l>l ^ (a;2-u;2)2+^4<^2 ^- ) 

of independent oscillators labeled by the index n. Each of these oscillators consists of a point mass 
w^n = [■0n(-f')]~^ attached to an anelastic spring with the complex spring constant kn{oS) = A;„[l + 
i(/)n(aj)], such that the resonant angular frequencies are tOn = \/h^rn^. So, in order to obtain the 
thermal noise spectrum one needs to find all the normal modes, their effective masses, resonant 
frequencies, and loss functions. 



4.3 Modes of a pendulum suspension 

The most important modes of a pendulum suspension are the pendulum mode, the rocking mode, and 

the violin modes. We will not consider here the rocking mode because for multi-loop suspensions the 
rocking motion of the test mass is essentially suppressed. The loss function of each mode depends on 
the type of mode and on anelastic properties of the pendulum wire. 



4.3.1 The pendulum mode 

For the pendulum mode, we will assume that the mass of the wire is much smaller than the mass of 

the bob (which is the test mass) and that the bob is attached near its center of mass. Also, the angle 
by which the pendulum swings is considered to be very small. Then the pendulum may be modelled 
as an oscillator of the resonant angular frequency 



a;p = .y^, (4.15) 



where g is the acceleration due to the Earth gravity field, and L is the pendulum length. 

The energy of the pendulum consists of two parts: the gravitational energy and the elastic 
energy £^el due to the bending of the wire. The gravitational energy is lossless; provided that all 
the losses due to interactions with the external world (friction in the residual gas, dumping by eddy 
currents, recoil losses into the seismic isolation system, friction in the suspension clamps, etc.) are made 
insignificant by careful experimental design, the assumption is made that the losses are dominated by 
internal friction in the wire material. Consequently, Af = A£"ei. Usually, £"gr ^ S^y^ so we obtain for 
the pendulum-mode loss function: 

(4.16) 



27r£^tot 27r(£^el + ^gr) 27r£^el £gr 

Note that 

is the loss function for the wire itself which occurs due to anelastic effects in the wire material. Then 
we obtain 

0P = Cp(pw, (4.18) 
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where .^p is the ratio between the elastic energy and the gravitational energy for the pendulum mode, 



The lossless gravitational energy of the pendulum is 

= ^MujlL^el = \MgLel, (4.20) 



where M is the pendulum mass and Om is the maximum angle of swing. The elastic energy depends 
on how many wires are used and how they are attached to the pendulum. For one wire, the fiber in 
the pendulum mode will bend mostly near the top, with the bending elastic energy 



^ei = \VTEiel. (4.21) 



Here, T is the tension force in the wire (T = Mg for one wire), E is the Young modulus of the wire 
material, and / is the moment of inertia of the wire cross section (/ = ^vrr^ for a cylindrical wire of 
radius r). Using these results, one finds for a single-wire pendulum: 




^/^^^^ (4.22) 

^ 2MgL nr \l nr \l rj. \ J 

This result can be easily generalized for the case when the test muss is suspended by wires. Then 



the elastic energy £c,\ of Eq. ( 4.21 ) should be multiplied by N and the tension in each wire becomes 
T = Mg/N. Then 

2MgL 2L\ Mg 2L^ T ^ ' 

In Eq. ( |4.23| ) we assumed that all the wires are in one plane: a plane through the center of mass of 
the pendulum, whose normal is parallel to the direction of swing. (Note that in such an configuration 
one should take into account the rocking mode of the test mass.) In this arrangement, the pendulum 
mode causes bending of the wires mostly at the top. If one uses a number of wire loops along the 
test mass length, then the rocking mode is essentially suppressed and the wires bend both at the top 
and the bottom. Therefore, the bending elastic energy of the multi-loop configuration is given by 



multiplying the result of Eq. (4.21) by 2N 



N 



fei = -jVTEiel. (4.24) 
Then the energy ratio is 

nVtEI 1 I BIN 1 [eI 

= -11^ = iV w = IV ^- ^^^^ 

The contribution of the pendulum mode to the thermal noise spectrum is obtained from Eq. ( [4.14| ) 
by taking m„ = M, /c„ = Mg/L, ujn = ojp and 0n = = ip4>w This gives 

= ^r^^%^- (4-26) 
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For LIGO suspensions, /p = io^jl-K is about 1 Hz. This is much below the working frequency range 
(near 100 Hz), so we may assume Wp/a; <C 1. Also, the loss function is very small, 0p < 10~^. Then 
the pendulum-mode contribution to the thermal noise spectrum is 



4.3.2 The violin modes 

The angular frequency of the nth violin mode (n = 1, 2, 3, . . .) is given by 

nvr Ft 



L V p 



2 1 / nvr \ 2' 
1 + ^ + 



(4.28) 



where L is the length of the wire, T is the tension force, p is the linear mass density of the wire, and 

ke - (4.29) 

In the violin mode the wire bends near both ends in a similar way. The bending occurs over the 
characteristic distance scale ~ yjEI /T, the same as in the pendulum mode. For ^ L, which 
is a very good estimation for heavily loaded thin wires like in LIGO, one have approximately, 

nvr It , , 

0..=.-^-. (4.30) 

This is just the angular frequency of the nth vibrational mode of an ideal spring. 
It can be shown that for the nth violin mode, the loss function is 

where the energy ratio is 

2 / n^vr^X 2 [eI ( 1 [eI , ,\ , 

Since keL ^ 1, for first several modes the energy ratio is approximately 

in^i. = \sl^. (4.33) 

This expression takes into account only the contribution to the elastic energy due to wire bending 
near the top and the bottom. For higher violin modes, one should also consider the contribution due 
to wire bending along its length, which leads to Eq. ( [4.32 ). 



For the one-loop suspension configuration, the elastic energy of the lowest violin modes is about 
twice of that for the pendulum mode (for the last one the wires bend only at the top while for the 
former ones the wires bend at both ends). In the multi-loop configuration, the elastic energy of the 
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lowest violin modes and of the pendulum mode is approximately the same. On the other hand, the 
gravitational energy of the pendulum mode is by a factor of 2 larger than that of a violin mode. For 
the violin modes of each wire, the gravitational energy is ^TL9'^. Then for N wires, 

(S^r). = \nTL01 = \MgLel. (4.34) 

This is just one half of the gravitational energy for the pendulum mode, (fgr)p = ^MgL9^ (cf. 
Eq. ( 4.20| )). This explains the difference between the loss functions for the pendulum mode and for 



the violin modes: — 4^p for the one-loop configuration and — 2^p for the multi-loop configuration. 
The effective mass of the nth violin mode is 

m„ = [ML)]-' = (^y^ (4.35) 

where we took expression (1^ for LOn and T = Mg/N. This effective mass arises because the violin 
vibrations of the wire cause only a tiny recoil of the test mass M. The contribution of the violin modes 
to the thermal noise spectrum is given by 

5;^.) = 1^£^=£!;^|^. (436) 



n= 



Typical values of /i = u;i/27r are from 350 to 500 Hz. If we are interested in the thermal spectrum 
density near 100 Hz, we can assume ^n- Then we have approximately 

n=l "■ n=l 

One can see that the contributions of higher violin modes are very small due to the factor in the 
sum. Taking 0„ = ^n4>w and using Eq. ( 4.32|) , we obtain 



E ^ = rr + ^w(^) - 7r7\hr^w(o.). (4.38) 



, n4 fce^\90 UkeL r ^ ' A5Ly T 

n=l \ / ' 



Here, we assumed k^L ^ 1. Finally, we substitute ( 4.38| ) into ( 4.37 ) and find the following expression 



for the violin-mode contribution to the thermal noise spectrum. 



(-) - T^^bTp L y (4.39) 

5 Experiments on anelasticity effects for pendulum suspensions 

5.1 Basic types of experiments 

In order to predict the thermal noise fluctuations in pendulum suspensions, two basic types of exper- 
iments are performed: 

1. Investigations of anelastic properties of wires made of various materials, in order to determine 
the wire loss function (^w('^)- 



13 



2. Measurements of quality factors {Q = 4> ^ at a resonance) for the pendulum and violin modes 
of actual suspensions, in order to verify the relationships 

0p(^) = ^p0w(w), (Av(w) = i^(t)^{uj). (5.1) 

Early experiments showed serious discrepancy between the measured quahty factors and those pre- 
dicted using Eq. ( |5.1D . It was discussed that this discrepancy may happen due to stress-dependent 
effects. However, it was shown later that the internal losses of the wires are almost independent of the 
applied stress. Many recent experiments proved that the above discrepancy appears due to serious 
losses in the clamps. A smart design of clamps can be used to reduce these excess losses and then 
predictions of Eq. ( ^.1] ) are quite accurate. A very promising possibility is the use of monolithic or 
semi-monolithic suspensions. The design of clamps plays a crucial role in the reduction of the thermal 
noise of the test mass suspensions. 

5.2 Internal losses in wire materials 

A number of experiments were performed to study internal losses of various wire materials (e.g., steel, 
tungsten, fused quartz, and some others). The main drawback of many of these experiments is a 
small number of frequencies for which was measured. Also, there are serious discrepancies between 
results of different experiments. Therefore, the exact frequency dependence of (j)^ is still unclear for 
many materials. Below, we briefly review results of some recent experiments. 

Kovalik and Saulson, 1993 Method: Quality factors were measured for resonances of freely sus- 
pended wires. Materials: Tungsten, silicon, sapphire, fused quartz. Results: Insignificant frequency 
dependence for tungsten; for fused quartz, measured cf)^ are above those predicted by the thermoelastic 
damping (TED) for some frequencies and near TED for others; sapphire and silicon showed behavior 
consistent with TED. 

Saulson et al., 1994 Method: Quality factors were measured for an inverted pendulum of tunable 
length. Material: Free-Flex cross-spring flexure made of crossed steel strips. Results: In agreement 
with a frequency-independent (p^. 

Gillespie and Raab, 1994 Method: Quality factors were measured for resonances of freely sus- 
pended wires. Material: Steel music wires. Results: A constant value of (p^ for low frequencies (from 
30 to 150 Hz). For higher frequencies (from 150 Hz to 2 kHz) increases with oj, like TED predicts, 
but the measured value (/>meas is well above 4>teb- These results may be explained by the formula 
^^meas = <^TED + 0ex5 where (/>ex is a frequency-independent excess loss. 

Rowan et al., 1997 Method: Quality factors were measured for resonances of ribbons fixed at one 
end. Material: Fused quartz ribbons. Results: Data were obtained for 5 resonances in the range from 
6 to 160 Hz. 0meas IS Well above i?i)TED for lower frequencies (below 30 Hz), and in agreement with 
TED for higher frequencies (above 80 Hz). 

Dawid and Kawamura, 1997 Method: Quality factors were measured for the violin modes of 
wires fixed at both ends in a "guitar"-type apparatus. Materials: Invar, titanium, steel, tungsten 
and several other metals. Results: (/"meas proportional to \/T, in accordance with the formula 
</)v = {2 / L) ^JEI/T 4>^ for frequency-independent (f)^. 
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Huang and Saulson, 1998 Method: Quahty factors were measured for resonances of freely sus- 
pended wires. Materials: Steel and tungsten. For steel, (pmeas coincides with the predictions of TED 
(the characteristic Debye-peak frequency dependence). Some differences were found between proper- 
ties of annealed wires (^meas slightly above ^ted) and "curly" wires (^meas slightly below ^ted)- The 
difference can be explained by modifications of thermal properties. For tungsten wires, c/'mcas only 
slightly increases with frequency; the loss function increases with the wire diameter, as should be for 
TED at frequencies well below /. 

6 Conclusions 

It is seen that predictions of the spectral density for thermal fluctuations in pendulum suspensions 
depend strongly on the type of the dissipation mechanism. Sources of external losses (friction in the 
residual gas, dumping by eddy currents, recoil losses into the seismic isolation system, friction in the 
suspension clamps, etc.) should be eliminated by careful experimental design. In particular, results 
of many recent experiments show that excess losses in clamps may seriously deteriorate the quality 
factors of suspension resonances. When external losses are made sufficiently small, the main source 
of dissipation is the internal friction in the wires due to anelastic effects. The thermal noise spectrum 
depends on the form of the loss function. Unfortunately, the exact frequency dependence of the wire 
loss function is not yet completely understood. In many experiments cp^ was measured only 

at few frequencies and experimental uncertainty of results was often quite large. Moreover, there is a 
contradiction between results of different experiments. Therefore, it is very difficult to make certain 
conclusions about the behavior of 4>^{uj). In particular, it is unclear if clamp losses are negligible in 
experiments with freely suspended wires, as is usually assumed. Certainly, there is a room for more 
experiments on anelastic properties of wires, in order to clarify the issue of internal friction in the 
frequency range of interest for gravitational- wave detection. 

Acknowledgment 

This work would be impossible without great help and encouragement by Malik Rakhmanov. I thank 
him for long hours of illuminating discussions and for encouraging me to enter the realm of thermal 
noise and anelasticity. 

Appendix: Correlation function and spectral density 

Consider a system characterized by some quantity a (e.g., position or velocity). For stationary pro- 
cesses, the correlation function is 

Pa{t) = (a(r)a(r + 0), 

where the average is over a statistical ensemble. Using the ergodic theorem, this can be replaced by 
the time average, 

1 

pa(t)= lim - / a{t')a{t' + t)dt' . 

T— >oo T J—T 

Now, define the function 




t G [-T,T] 
other 
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and its Fourier transform 

/'OO 
-oo 

The definition of the spectral density is 

|aT(w)P 
Sa{u;) = lim ■ 

T— >oo TTi 

It is easy to see that the correlation function Pa{t) and the spectral density Sa{io) are related via the 
Fourier transform: 

1 fOO 1 poo 

Pa{t) = / SMe-'^'(kj, SaH = - / Pa{t)e"^'dt. 

^ J —oo J —oo 

This result is known as the Wiener-Khinchin theorem. 
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